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Abstract 

A  two-dimensional,  steady  state  model  for  proton  exchange  membrane  fuel  cell  (PEMFC)  is  presented.  The  model  is  used  to  describe  the 
effect  of  flow  mode  (coflow  and  counterflow),  operation  conditions  and  membrane  thickness  on  the  water  transport,  ohmic  resistance  and 
water  distribution  in  the  membrane,  current  density  distribution  along  the  channel  and  performance  of  PEMFC.  Effect  of  liquid  water  on 
the  transport  in  the  two-phase  region  of  cathode  diffusion  layer  was  considered.  Water  transport  in  the  membrane  by  electro-osmosis  drag, 
diffusion  and  convection  were  combined  in  this  model.  The  model  predicts  that  the  dry  reactant  gases  can  be  well  internally  humidified  and 
maintain  high  performance  when  PEMFC  is  operated  in  the  counterflow  mode  without  external  humidification.  Counterflow  mode  does 
not  show  any  advantageous  while  the  reactant  gases  are  high  humidified  or  saturated.  Compared  to  the  coflow  mode,  counterflow  mode 
improves  the  current  density  distribution  with  dry  or  low  humidity  gases.  The  higher  the  anode  is  humidified,  the  more  water  will  migrate 
from  anode  to  cathode.  The  modeling  results  compare  very  well  with  experimental  results. 
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1.  Introduction 

Proton  exchange  membrane  fuel  cells  (PEMFCs)  are  ef¬ 
ficient  and  environmentally  clean  electrical  generators  that 
are  being  developed  for  both  stationary  and  mobile  applica¬ 
tions  [1-4].  Over  the  past  decade,  significant  research  effort 
of  PEMFC  has  been  focused  on  fuel  cell  structure  design, 
the  development  of  better  catalyst,  the  better  structure  of 
catalyst  layer  and  gas  diffusion  layer,  the  development  of 
membrane  with  high  performance,  thermal  and  water  man¬ 
agement  [5]. 

In  the  PEMFC  application,  high  membrane  conductivity 
is  one  of  the  key  parameters  for  attaining  high  performance. 
The  most  commonly  used  exchange  membrane  in  fuel  cell  is 
Nafion  (DuPont).  The  proton  conductivity  of  Nafion  mem¬ 
brane  increases  with  the  water  content  [6,7].  In  order  to 
maintain  high  water  content  of  membrane  along  the  channel, 
the  reactant  gases  are  humidified  by  passing  them  through 
humidifiers.  High  humidity  of  reactant  gases  may  lead  to  liq¬ 
uid  water  existed  in  the  cathode.  When  cathode  is  flooded, 
the  diffusion  of  reactant  gas  through  its  pores  is  difficult. 
For  the  H2/O2  system,  water  flooding  is  known  to  be  the 
main  limiting  factor  that  restricts  the  improvement  of  per¬ 
formance  at  high  current  densities  [8].  In  addition,  external 
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humidification  of  the  reactant  gases  is  a  burden  for  the  fuel 
cell  system.  In  order  to  simplify  the  subsystem  and  reduce 
heat  removal  from  the  PEMFC  stacks,  operation  of  PEM¬ 
FCs  without  external  humidification  of  the  reactant  gases 
is  advantageous  for  the  PEMFC  system  [9].  Qi  designed  a 
double-path-type  flow-field  that  has  two  gas  inlets  and  two 
gas  outlets  for  PEMFCs.  Such  a  design  enables  the  entering 
dry  reactant  gases  to  be  dydrated  by  the  exiting  gas  [10]. 
The  problem  of  water  transport  in  PEMFC  with  different 
humidification  conditions  has  been  the  object  of  several  ex¬ 
perimental  and  theoretical  studies  [6,11-15]. 

Models  for  water  management  are  useful  for  describing 
and  getting  great  insight  into  the  mass  transport  processes. 
One-dimensional  PEMFC  models  were  developed  in  the 
early  1990s  by  Bernardi  and  Verbrugge  [14]  and  Springer 
et  al.  [13].  In  the  hydraulic  model  developed  by  Bernardi 
and  Verbrugge  [14],  the  author  assumes  that  the  membrane 
is  fully  saturated  with  water  and  that  most  of  the  water  is 
transported  through  the  electrodes  in  the  liquid  phase.  This 
model  may  be  not  applicable  in  the  case  of  feeding  dry  or 
unsaturated  reactant  gases.  Springer  et  al.  [13]  presented 
a  diffusion  model.  Water  transport  in  the  membrane  by 
means  of  electro-osmotic  drag  and  back  diffusion  process  at 
different  humidification  conditions  was  described  in  detail 
in  the  diffusion  model.  Springer  et  al.  [6]  assumed  that  the 
presence  of  liquid  water  only  modifies  the  effective  porosity 
for  oxygen  gas  transport  whereas  Bernardi  assumed  that  a 
constant  liquid  volume  fraction  and  no  interactions  between 
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Nomenclature 

a  electrode  length  (m) 

b  electrode  width  (m) 

c  concentration  of  water  in  the  membrane 

(mol/m3) 

Cf  fixed  charge  concentration  in  membrane 
(mol/m3) 

D  diffusion  coefficient  (m2/s) 

Dm  diffusion  coefficient  of  water  in  the  membrane 
(m2/s) 

Eq  open  circuit  potential  (V) 

£ceii  cell  potential  (V) 

/  local  fraction  of  the  water  produced 

that  evaporates  through  the  anode 
F  Faraday  constant  (96,485  C/mol) 

I  local  current  density  (A/m2) 

Io  exchange  current  density  (A/m2) 

kp  hydraulic  permeability  of  water  in 

membrane  (m2/s) 
m  molecular  weight  (g/mol) 

M  mole  flow  (mol/s) 

rid  electro-osmotic  drag  coefficient 

N  mole  flux  (mol/(m2  s)) 

p  pressure  (Pa) 

Pbm  logarithmic  average  value  of  nitrogen 

partial  pressure  (Pa) 

R  universal  gas  constant  (8.314  J/(molK)) 

Rco  contact  resistance  (Q  m2) 

Rm  membrane  resistance  (Q  m2) 

.v  liquid  water  saturation 

T  temperature  (K) 

v  direction  along  the  channel  length  (m) 

Xi  mole  fraction  of  species  i 

y  direction  across  the  cell  (m) 

Greek  letters 

a  water  relative  humidity 

8  thickness  (m) 

8ij  Angstrom  unit 

dimensionless  function  of  the  temperature  and 
of  the  intermolecular  potential  field  for  the 
molecules  of  i  and  j 
s  porosity 

r]  electrode  overpotential  (V) 

X  water  content  of  membrane 

fi  water  viscosity  (Pa  s) 

a  membrane  conductivity  (£2m)-1 

§  stoichiometric  flow  ratio  of  reaction  gas 

Superscripts  and  subscripts 
0  initial  value 

a  anode 

av  average 


c  cathode 

d  diffusion 

e  electro-osmotic  drag 

H  hydrogen 

ij  species  i,j 

m  membrane 

O  oxygen 

p  pressure 

w  water 


liquid  and  gas  flow.  One-dimensional  model  is  so  simple 
that  it  neglects  the  variations  of  water  content  of  membrane 
and  gas  concentration  along  the  flow  channel.  It  cannot 
be  used  to  optimize  operation  conditions  and  to  predict 
performance  of  PEMFC  with  high  working  surface  area. 

To  investigate  and  optimize  the  distribution  of  gas  com¬ 
ponent,  water  content  of  membrane,  current  density  and 
net  water  transport  along  the  flow  channel,  two-dimensional 
models  were  developed  in  recent  years  [15-23].  Fuller  and 
Newman  [15]  simulated  a  PEMFC  operation  under  the  con¬ 
dition  of  moist  gas  in  electrodes.  In  this  model,  concentrated 
solution  theory  is  used  for  description  of  the  membrane. 
The  water  distribution  in  the  membrane  and  the  net  water 
flux  across  it  were  obtained.  Nguyen  and  White  [16]  com¬ 
pared  three  humidification  designs  using  a  combined  heat 
and  mass  transfer  model.  Both  Fuller  and  Nguyen  consider 
coflow  pattern  to  study  transfer  processes. 

Recently,  Gurau  et  al.  [17]  presented  a  comprehensive 
model  for  the  entire  sandwich  of  a  PEMFC  including  the 
gas  channels.  Similarly,  Hsing  and  Futerko  [18]  developed 
a  two-dimensional  model  of  coupled  fluid  flow,  mass  trans¬ 
port  and  electrochemistry  of  PEMFC  operating  without  ex¬ 
ternal  humidification  of  hydrogen  and  air  streams.  Hsing 
assumed  a  constant  local  current  density  in  the  model.  Dan- 
nenberg  et  al.  [19]  presented  an  along-the-channel  mass  and 
heat  transfer  model  for  a  PEMFC  and  an  agglomerate  ap¬ 
proach  was  used  for  the  description  of  the  active  layer  of 
the  cathode.  Yi  and  Nguyen  [20]  analyzed  multicomponent 
transport  in  porous  electrodes  of  PEMFCs  using  the  inter- 
digitated  gas  distributors.  It  should  be  noted  that  all  these 
models  do  not  account  for  the  mass  transport  processes  in 
the  present  of  liquid  water. 

During  PEMFC  operation,  especially  at  low  stoichiomet¬ 
ric  flow  rate  of  reactant  gases,  liquid  water  is  likely  to  appear 
in  the  cathode,  resulting  in  two-phase  transport  phenomena. 
Wang  et  al.  [21]  modeled  the  two-phase  flow  and  transport 
in  the  air  cathode  of  PEMFCs  and  concluded  that  capillary 
action  was  the  dominant  mechanism  for  water  transport 
inside  the  two-phase  zone  of  the  hydrophilic  structure.  The 
transport  processes  become  significantly  complicated  in  the 
two-phase  zone.  A  two-dimensional,  two-phase,  multicom¬ 
ponent,  transient  model  was  developed  by  Natarajan  for  the 
cathode  of  the  PEMFC  [22].  It  was  found  that  the  perfor¬ 
mance  of  the  cathode  was  dominated  by  the  dynamics  of 
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liquid  water,  especially  in  the  high  current  density  range. 
Neglecting  the  transport  of  dissolved  reactant  gases,  this 
liquid  phase  in  the  electrode  is  a  one-component  system. 
Janssen  [23]  reported  that  the  presence  of  a  gradient  in  the 
chemical  potential  of  liquid  water  could  produce  a  flux. 

The  aim  of  the  present  work  was  to  develop  a  model  that 
describes,  in  steady  state,  the  water  transport  in  the  mem¬ 
brane  and  the  fuel  cell  performance,  at  dry  and  different  hu¬ 
midified  gases  in  the  coflow  and  counterflow  modes.  Effect 
of  liquid  water  on  the  effective  porosity  for  gas  transport  in 
the  cathode  was  considered  in  this  model.  Such  a  study  will 
be  useful  for  the  design  and  operation  of  practical  PEMFC 
stacks. 
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Fig.  1.  A  schematic  of  PEMFC  regions:  (1)  anode  flow  channel;  (2)  gas 
diffusion  anode;  (3)  membrane;  (4)  gas  diffusion  cathode;  (5)  cathode 
flow  channel. 


2.  Experimental  and  model  development 

Membrane  electrode  assemblies  (ME A)  of  140  cm1 2  were 
prepared  using  Nation  112  and  Nation  117  membranes 
(DuPont)  and  electrodes  made  by  our  groups  containing 
0.4  mg/cm2  platinum  (20%  Pt/Vulcan  XC-72  carbon  pow¬ 
der,  Cabot  Corp.)  The  membrane  was  first  pretreated  using  a 
usual  procedure,  i.e.  maintained  in  3%  H2O2  at  353  K  for  1  h, 
then  washed  with  pure  water.  The  membrane  was  then  main¬ 
tained  in  0.5  mol/l  H2SO4  at  353  K  for  1  h  and  washed  with 
pure  water.  Electrodes  were  impregnated  with  recast  Nation 
solution  (Aldrich)  by  brushing.  The  ME  A  was  prepared  by 
hot-pressing  a  membrane  sheet  between  the  two  electrodes. 

Pure  hydrogen  and  oxygen  were  used  as  reactant  gases, 
and  stoichiometry  of  hydrogen  and  oxygen  were  kept  at 
needed  values.  The  electrical  parameters  of  the  cell  were 
controlled  by  an  automatically  test  station  interfaced  to  a 
personal  computer.  The  amount  of  water  exiting  the  cell  at 
the  anode  and  cathode  gas  outlets  was  measured  by  con¬ 
densing  the  water  in  an  ice  trap  and  then  by  weighting  the 
condensed  water.  For  water  transport  measurements,  the  cell 
was  operated  in  the  constant-current  mode. 

2.7.  Model  assumptions 

A  two-dimensional,  steady  state,  water  transport  model 
for  a  PEMFC  is  developed.  A  schematic  illustration  of  a 
PEMFC  in  the  counterflow  mode  is  shown  in  Fig.  E  The 
model  regions  consist  of  a  membrane  sandwiched  between 
two  gas  diffusion  electrodes,  the  anode  flow  channel  and 
cathode  flow  channel.  The  model  accounts  for  the  water 
content  of  membrane,  the  mass  transport  of  water  and  other 
gaseous  species  across  the  cell  (y-direction),  and  current 
density  distribution  along  the  flow  channel  (x-direction,  di¬ 
rection  of  the  hydrogen  flow).  To  keep  the  water  transport 
model  simple,  a  few  assumptions  are  made. 

1 .  The  temperature  is  assumed  to  be  uniform  throughout  the 

membrane  electrode  assembly  and  gas  channels. 

2.  An  ideal  gas  mixture  is  assumed.  In  the  flow  channel, 

a  plug-flow  condition  is  assumed,  and  pressure  drop 


along  the  channel  is  neglected.  The  gases  and  liquid 
water  are  treated  as  being  perfectly  well  mixed  in  the 
y-direction. 

3.  The  diffusion  layer  is  uniform  and  the  mass  transport  in 
the  electrodes  is  considered  only  along  the  y-direction. 

4.  It  is  assumed  that  the  catalyst  layer  is  ultra-thin,  gas  dif¬ 
fusion  through  the  catalyst  layer  is  assumed  to  be  negli¬ 
gible. 

5.  Water  contents  of  membrane  on  the  interface  of  both  sides 
are  assumed  to  be  in  equilibrium  with  water  vapor  in 
gases.  The  mass  transport  of  hydrogen  and  oxygen  in  the 
membrane  is  negligible.  Water  transport  in  the  membrane 
is  considered  only  along  the  y-direction. 

6.  It  is  assumed  that  liquid  water  transport  is  not  a  limited 
process.  How  liquid  water  is  transferred  in  the  cathode  is 
not  included  in  this  model.  It  is  assumed  that  the  change 
in  gas  pressure  caused  by  capillary  effects  is  small  enough 
to  have  a  negligible  effect  on  gas  transport. 

7.  Electrode  pores  for  gas  flow  are  separated  from  pores  for 
liquid  water.  Liquid  water  modifies  the  effective  porosity 
for  oxygen  gas  transport  in  the  cathode  diffusion  layer 
and  affects  the  cathode  overpotential. 

2.2.  Mass  balance 

There  are  three  water  transport  mechanisms  in  the  mem¬ 
brane  [24] :  (i)  electro-osmotic  drag,  which  is  caused  by  pro¬ 
ton  transport  from  the  anode  side  to  the  cathode  side  of  the 
membrane;  (ii)  back-diffusion  by  the  concentration  gradient 
of  water,  the  concentration  gradient  is  due  to  the  transport 
of  water  by  the  electro-osmotic  drag  and  the  production  of 
water  at  the  cathode  via  the  oxygen  reduction  reaction;  and 
(iii)  convection  generated  by  pressure  gradient  within  the 
pores  of  the  proton  exchange  membrane. 

In  a  calculated  unit  of  membrane,  there  is  no  potential 
gradient  in  the  x-direction.  The  electro-osmotic  drag  flux  in 
y-direction,  is  as  follows 

7(x) 

Wv  ,y,Q  —  nd 


F 


(1) 
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where  n&  is  the  electro-osmotic  drag  coefficient  (number  of 
water  molecules  carried  by  a  proton),  7(x)  the  local  current 
density  of  the  fuel  cell,  and  F  the  Faraday’s  constant.  There  is 
no  pressure  gradient  within  the  membrane  in  the  x-direction, 
so  convection  generated  by  pressure  gradient  is 


p 


kp  dp 
— -c  — 
M  dy 


where  kv,  /jl,  p  and  c  are  the  hydraulic  permeability  of  water 
in  membrane,  water  viscosity,  pressure  and  concentration 
of  water  in  the  membrane,  respectively.  Diffusion  flux  is  as 
follows 


Nw,y,d 


3  c 


N W,A',cl  =  —D 


m 


3  c 
dx 


where  Dm  is  the  diffusion  coefficient  of  water.  As  the  thick¬ 
ness  of  membrane  is  much  less  than  its  length  (y  is  less  than 
x  by  several  orders),  the  diffusion  of  water  in  the  membrane 
in  x-direction  can  be  neglected.  The  calculated  results  will 
show  that  concentration  gradient  of  water  in  y-direction  is 
higher  several  orders  than  that  in  x-direction.  n d  and  Dm 
were  measured  by  Springer  et  al.  [13] 


n  d 


2.5A 

~72 


Dm  = 


10  10  exp 


2416 


_LWI 

303 )  \T 

2 


x  (2.563  -  0.33A  +  0.0264AZ -0.00067^),  A  >  4 

(6a) 

Eq.  (6a)  is  the  fit  for  A  >  4,  as  for  2  >  A  <  4,  Dm  can  be 
expressed  as 


Dm  ~  10  10  exp 


2416 


1 


303 


1 

T 


x  (— 1.25A  +  6.65),  3  >  A  <  4 


(6b) 


I)m  —  10  10  exp 


x  (2.05A  -  3.25),  2  >  A  <  3 


where  A  is  the  water  content  of  membrane  (water  molecule 
per  SOzT).  According  to  Bernardi  and  Verbrugge  [25],  the 
pressure  profile  through  the  membrane  region  is  linear,  and 
the  velocity  is  constant. 

dp  =  Pc  ~  Pa  m 

dy  ~  Sm  { 

Water  concentration  is  calculated  as  follows 


where  Cf  is  the  fixed  charge  concentration  in  membrane. 
The  water  flux  through  the  membrane  (y-direction)  can  be 
written  as  follows 


2.5A  f(x) 
~~22  F7” 


In  the  anode,  gas  diffusion  in  x-direction  is  much  less  than 
that  in  y-direction,  and  therefore  gas  diffusion  in  x-direction 
is  neglected.  Stefan-Max  well  equation  is  used  to  describe 
anode  mass  transport. 


d  Xj 

dy 


(10) 


where  DfJf  is  an  effective  binary  diffusivity  of  the  pair  i-j,  T 
the  temperature,  and  R  the  universal  gas  constant.  We  relate 
the  effective  diffusion  coefficient  to  the  diffusion  coefficient 
in  a  nonporous  system  by  [25] 


Df  =  Djjs15 


(11) 


Binary  diffusion  coefficients  can  be  calculated  by  us¬ 
ing  the  equation  commended  by  Bird  et  al.  [26] 


Dij  =  1.8583  x  1CT7 


yr3[(l/m,0  +  (1/m;)] 
(p/po)Sjj^D 


(12) 


where  m*,  mj  are  the  molecular  masses  of  i  and  j,  respec¬ 
tively,  8ij  the  Angstrom  unit,  ^2d  the  dimensionless  function 
of  the  temperature  and  of  the  intermolecular  potential  field 
for  the  molecules  of  i  and  j. 

At  steady  state,  the  fluxes  of  all  the  diffusing  species  in 
the  pores  are  constant,  and  the  amount  of  water  produced  by 
the  electrochemical  reaction  is  equal  to  the  amount  of  water 
transported  through  the  pores  of  diffusion  layer.  The  flux  of 
hydrogen,  oxygen  and  water  can  be  described  [27] . 

M.  =  —2Nq  =  Nh  =  -JV*  +  Nl  (13) 


The  parameter  /  is  defined  as  the  local  fraction  of  the 
water  produced  that  evaporates  through  the  anode  layer. 


(14) 


l  ra 

/a V  =  —  /  fldx  (15) 

v  Jo 


where  a  is  the  length  of  channel,  7av  the  average  current 
density.  For  two  components  (H2O,  H2),  Eq.  (10)  can  be 
solved.  The  concentration  distribution  of  water  in  the  anode 
diffusion  layer,  x^(y),  and  anode  hydrogen  flux  are  obtained 
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Eqs.  (16)  and  (17)  are  not  fit  for  /  =  1.  As  for  /  =  1,  x^(y) 
and  Ah  are 

=  |-(-C,a  _xw)  +  xw>  /  =  i  (18) 

D  Z)eff 

^  =  ^(<•“-4).  /  =  1  (19) 

At  the  inlet  region  of  cathode,  humidified  or  dry  oxygen 
enters  the  fuel  cell,  and  no  liquid  water  exits  in  the  flow 
channel  and  electrode.  Stefan-Maxwell  equation,  Eq.  (8),  is 
used  to  describe  cathode  mass  transport  in  this  single-phase 
region. 

While  oxygen  flow  along  the  flow  channel,  oxygen  is 
humidified  by  water  produced  by  electrochemical  reaction. 
Liquid  water  exits  in  the  cathode  as  oxygen  is  saturated  and 
effect  of  liquid  water  on  the  transport  in  the  two-phase  region 
of  cathode  diffusion  layer  is  considered.  Electrode  pores  for 
gas  flow  are  assumed  to  separate  from  pores  for  liquid  water. 
In  the  pore  of  diffusion  layer  for  gas  flow,  steam  pressure  is 
saturated  and  concentration  gradient  of  inert  impurity  gases, 
e.g.  nitrogen,  is  formed.  It  is  obvious  that  the  nitrogen  flux 
is  zero  in  steady  state.  The  oxygen  flux  though  the  cathode 
diffusion  layer  in  two-phase  region  is  described  as  follows: 


RT8cpbm 


Po)  and  PBM 


P  B2  —  /’Bl 
ln(/?B2/PBl) 
(20) 


where  §  is  the  stoichiometric  flow  ratio  of  reaction  gas,  the 
ratio  of  the  amount  of  reactant  gases  added  to  the  fuel  cell  to 
the  amount  that  is  required  by  the  electrochemical  reaction;  b 
is  width  of  electrode.  On  the  counterflow  operation,  oxygen 
and  water  mole  flow  rates  are 


Mo(x)  =  (£c  -  1)*  [  Nodx  +  bf  No  dx, 

Jo  Jo 

(counterflow)  (26) 


Mcw(x)  =  Mcw0  +  b[  (1  —  f)Nu  dx 

Jo 


'X 


—  b  I  (1  —  f)Nudx ,  (counterflow) 

Jo 


(27) 


A  fit  of  the  relationship  of  A  versus  water  relative  humidity 
(water  vapor  activity),  a ,  used  in  the  model  is  [13] 

X  =  0.043  +  17. 81a  —  39.85a2  +  36.0a3,  a  <  1  (28a) 


X  =  14.0  +  1.4(a  —  1),  1  <  a  <  3 

X  =  16.8,  a  >  3 
P 


a  =  x 


w 


p 


sat 


(28b) 

(28c) 

(29) 


2.3.  Cell  potential  and  current  density 


The  cell  voltage  (Eceii)  is  calculated  from  the  resis¬ 
tance  (include  membrane  resistance  and  contact  resis¬ 
tance,  Rco ),  current  density  and  electrode  overpotentials,  as 
follows: 

£ceii  =  Eo  -  rj(x)  -  ( Rm  +  Rco)I(x)  (30) 


Pbm  is  the  logarithmic  average  value  of  nitrogen  partial  pres¬ 
sure  on  both  side  of  cathode  diffusion  layer.  Liquid  water 
modifies  the  effective  porosity  for  oxygen  gas  transport  in 
the  cathode  diffusion  layer  [22]. 

£>^W  =  D0-W[fi(l-^)]1-5  (21) 


where  s  is  the  liquid  water  saturation.  The  presence  of  liquid 
water  in  the  porous  gas  diffusion  layer  results  in  mass  trans¬ 
port  limitations  imposed  by  the  restriction  of  gas  diffusion 
pores  [22]. 

While  reactant  gases  flow  along  the  flow  channel  (gas 
chamber  of  fuel  cell),  hydrogen  or  oxygen  transport  from 
channel  to  diffusion  layer  and  water  transport  from  diffusion 
layer  to  channel.  In  the  coflow  mode,  it  can  be  obtained  that 


ea  ex 

Mu(x)  =  §a b  /  Ah  dx  —  b  Ah  dx 

Jo  Jo 

eci  ex 

Mq(x)  =  £c b  Nodx  -  b  N0  dx, 

Jo  Jo 


M*(x)  =  M^0  +  bl  JNudx 
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(coflow)  (23) 
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Jo 
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In  the  most  of  model,  the  contact  electrical  losses  at  the 
interfaces  between  different  fuel  cell  elements  are  neglected 
[6,15-17].  Contact  resistance  is  taken  into  account  in  this 
model. 

A  fit  of  the  relationship  of  membrane  conductivity  versus 
X  is  [13] 


crT  —  exp 


1268 


1 


303 


x  (0.5139A  —  0.326),  for  A  >  1 


(31) 


While  X  is  less  than  1,  the  membrane  conductivity  was  as¬ 
sumed  constant.  The  membrane  resistance  is 


Rm  — 


-f 


’m 


1 


(T  t 


dy 


(32) 


Anode  overpotential,  when  pure  hydrogen  is  used,  is  small 
and  is  neglected  in  the  calculation.  Electrode  overpotential 
in  the  model  only  includes  cathode  overpotential.  77 (x)  can 
be  calculated  from  the  following  equation  [21]. 


RT 

05 ~F 
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I 

Jo(pc0/po)(l 


(33) 


here  To  is  the  exchange  current  density.  The  model  exchange 
current  density  parameter  was  adjusted  to  yield  results 
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that  suitably  mimic  the  experimental  results.  Liquid  water 
formed  at  the  interface  can  cover  the  active  reaction  sites. 
It  is  assumed  that  the  active  reaction  sites  decrease  linearly 
with  increase  in  the  liquid  water  saturation.  To  account  for 
this  coverage  phenomenon  the  (1  —  s)  term  is  include  in 
Eq.  (33).  The  average  current  density  for  a  given  cell  voltage 
is  calculated  by  integrating  along  the  whole  channel  length. 

1  fa 

4v  =  -  /  Kx)  dx  (34) 

a  Jo 

2.4.  Solution  technique 

The  inlets  of  hydrogen  and  oxygen  are  in  the  same  site  in 
coflow  mode.  The  model  equations  are  solved  numerically 
at  a  given  value  of  average  current  density,  7av.  Based  on  7av, 
mole  flow  rate  of  hydrogen,  oxygen  and  water  are  calculated 
in  the  inlet  of  gas  channel.  The  ECeii  and  a  values  of  I(x) 
are  initially  guessed,  and  one-dimensional  equations  for  the 
mass  transport  in  electrodes  and  membrane  are  solved.  Then 
/  and  distribution  of  water  content  are  obtained.  The  local 
current  density  is  calculated  and  compared  with  the  specified 
value.  If  the  calculated  value  disagrees  with  the  specified 
one,  a  new  value  of  I(x)  is  given.  This  process  is  repeated 
until  the  calculated  I(x)  agrees  with  the  specified  one.  This 
procedure  is  then  repeated  at  v  =  dv,  and  until  in  the  whole 
flow  channel.  Therefore,  7av  is  calculated  using  Eq.  (34)  and 
compared  with  the  initial  given  value.  A  value  of  Eceu  is 
given  again  and  secant  method  is  used  to  help  search  for  the 
correct  £cen  value. 

The  solution  process  with  counterflow  mode  is  similar  to 
that  with  coflow  mode.  The  site  of  the  inlet  of  hydrogen 
is  the  outlet  of  oxygen.  A  mole  flow  rate  of  water  in  the 
cathode  at  v  =  0  is  assumed.  The  mole  flow  rate  of  water  in 
cathode  at  v  =  a  is  calculated  and  compared  with  the  actual 
value.  The  mole  flow  rate  of  water  in  cathode  at  v  =  0  is 
obtained  using  secant  method. 

3.  Results  and  discussion 

Correct  water  management  has  been  found  to  be  impor¬ 
tant  to  maintain  high  proton  conductivity  of  membrane  and 
achieve  high  performance  for  PEMFCs.  A  mathematical 
model  of  PEMFC  is  very  useful  and  effective  tool  to  opti¬ 
mize  cell  designs  and  operation  conditions.  In  this  section, 
this  model  will  be  used  to  evaluate  flow  modes  and  vari¬ 
ous  operation  conditions  (relative  humidity,  temperature  and 
pressure  differential)  and  to  investigate  the  effect  of  mem¬ 
brane  thickness  on  the  performance  of  a  PEMFC.  Table  1 
lists  the  important  parameters  of  electrode  and  membrane 
used  as  the  base  case. 

3.1.  Comparison  of  the  calculated  and  experimental  results 

Fig.  2  compares  the  calculated  fuel  cell  potential  as  a 
function  of  current  density  with  experimental  data.  Electrode 


Table  1 

Parameter  values  of  electrode  and  membrane  used  as  the  base  case 


Parameter 

Value 

-^cell 

1.05  V 

7o 

40  A/m2 

Cf 

1.2  x  103  mol/m3* 

kp 

1.8  x  10-18  m2* 

[i 

3.56  x  10“6PaS* 

5m 

Nation  117 

1.75  x  10-4  m 

Nation  115 

1.25  x  10-4  m 

Nation  1135 

0.80  x  10“4m 

Nation  112 

0.50  x  10“4m 

a 

0.12m 

K 

3.0  x  10-4  m 

£a 

0.4 

5c 

3.0  x  10-4  m 

£c 

0.4 

Rco 

1  x  10-5  Qm2 

*  [14]. 


area  of  140  cm2  and  Nation  112  membrane  are  used  in  this 
test.  The  fuel  cell  is  operated  with  dry  gases  in  the  counter¬ 
flow  mode.  It  can  be  found  that  modeling  results  compare 
very  well  with  the  experimental  results. 

We  have  measured  the  fraction  of  product  water  leaving 
the  anode  side  of  a  PEMFC  over  a  wide  range  of  operation 
conditions.  Table  2  shows  the  comparison  of  the  calculated 
and  experimental  values  for  /av  at  different  operation  condi¬ 
tions.  Nation  117  membrane  is  used  in  this  test.  The  fuel  cell 
is  operated  in  coflow  mode.  The  temperature  of  fuel  cell  is 
kept  constant  at  353  K.  It  can  be  found  that  calculated  val¬ 
ues  for  /av  agreed  with  the  experimental  results.  Tests  1-3 
show  that  water  transfer  from  anode  to  cathode  increases 
with  increasing  relative  humidity  of  feeding  hydrogen. 

3.2.  Effect  of  relative  humidity  of  reactant  gases 

In  the  following  sections,  the  effect  of  various  cell  oper¬ 
ating  parameters  and  membrane  thickness  are  investigated 


Fig.  2.  Comparison  between  the  model  predictions  and  experimen¬ 
tal  results  when  Nation  112  membrane  is  used  (T  =  353  K, 
Pa  =  Pc  =  3  x  105  Pa,  £a  -  1-25,  £c  -  2.0,  aa  =  ac  =  0). 
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Table  2 


Calculated  and  experimental  values  for  /av 


Test 

fa 

fc 

Oia 

ac 

pa  (x  105  Pa) 

pc  (x  105  Pa) 

lay  (A/cm2) 

/av,exp 

fay, cal 

1 

1.20 

1.95 

0.61 

0.50 

4.0 

6.0 

0.402 

-0.078 

-0.072 

2 

1.20 

1.50 

1.01 

0.50 

4.0 

6.0 

0.410 

-0.146 

-0.141 

3 

1.20 

1.50 

1.63 

0.50 

4.0 

6.0 

0.403 

-0.272 

-0.265 

4 

1.20 

1.50 

1.06 

0.50 

4.0 

6.0 

0.606 

-0.160 

-0.157 

5 

4.01 

1.97 

0.20 

0.50 

4.0 

4.0 

0.175 

0.305 

0.241 

using  this  model.  Table  3  lists  the  cell  operating  parameters 
for  the  base  case.  Nafion  112  membrane  is  chosen  unless  it 
is  pointed  out  otherwise  in  the  text.  In  the  simulations,  one 
parameter  is  varied  and  the  other  parameters  are  kept  con¬ 
stant. 

Fig.  3(a)  and  (b)  shows  the  cell  potential  as  a  function  of 
current  density  for  various  relative  humidities  of  inlet  gases 
in  the  coflow  and  counterflow  modes,  respectively.  Cell  per¬ 
formance  increases  as  the  relative  humidity  of  inlet  reactant 
gases  increases  from  0  to  75%,  due  to  the  higher  water  con¬ 
tent  of  membrane  at  higher  relative  humidity  of  inlet  gases 
and  lower  membrane  resistance.  At  average  current  den¬ 
sity  7av  =  1.0  A/cm2,  the  cell  potential  varies  from  0.709  V 
(aa  =  ac  =  0.25)  to  0.723  V  (for  aa  =  ac  =  0.75)  in  the 
coflow  mode,  and  the  cell  potential  varies  from  0.698  V  (dry 
gas,  for  aa  =  ac  =  0.0)  to  0.723  V  (for  aa  =  ac  =  0.75)  in 
the  counterflow  mode.  The  model  predicts  that  the  cell  can¬ 
not  be  operated  with  dry  gases  in  the  coflow  mode.  While 
the  relative  humidity  of  inlet  reactant  gases  is  higher  than 
75%,  the  performances  of  fuel  cell  are  almost  constant  both 
in  the  coflow  and  counterflow  modes.  The  reason  why  cell 
performances,  at  high  humidity  of  inlet  gases,  are  rather  sim¬ 
ilar  is  that  the  decrease  of  membrane  resistance  is  compen¬ 
sated  by  increase  of  cathode  overpotential  due  to  the  liquid 
water.  To  reduce  heat  removal  and  liquid  water  in  the  cath¬ 
ode,  relative  humidity  of  inlet  reactant  gases  should  not  be 
higher  than  100%.  It  can  be  found  that  counterflow  mode 
cannot  show  any  advantageous  while  the  reactant  gases  are 
high  humidified  or  saturated. 

Fig.  4  shows  the  water  content  of  membrane  in  the  coun¬ 
terflow  mode.  Current  density  is  kept  constant  at  1.0  A/cm2 
and  the  fuel  cell  is  operated  with  dry  gases.  The  simulation 
results  show  that  water  content  gradient  across  the  mem¬ 
brane  is  high  in  the  inlet  region  of  hydrogen  and  that  the 

Table  3 


Operating  parameter  for  the  base  case 


Parameter 

Value 

la\ 

1.0  A/cm2 

Pa 

4  x  105  Pa 

Pc 

4  x  105  Pa 

T 

353  K 

(X  a 

0.5 

ac 

0.5 

fa 

1.5 

fc 

2.0 

water  content  of  membrane  is  very  low  in  the  outlet  region 
of  hydrogen.  It  can  be  obtained  that  the  water  content  gradi¬ 
ent  in  y-direction  is  2  orders  higher  than  that  in  v-direction. 
Therefore,  it  is  correct  to  neglect  the  diffusion  of  water  in 
the  membrane  in  y-direction. 

Fig.  5  shows  the  distribution  of  current  density  along  the 
gas  channel  for  various  relative  humidities  of  inlet  gases. 
Cell  potential  is  keep  constant  at  0.7  V  and  average  current 
density  may  be  varied.  In  the  inlet  region  of  gases,  the  cur¬ 
rent  density  is  very  low  with  low  humidity  reactant  gases  in 
the  coflow  mode.  It  is  due  to  the  low  water  content  of  mem¬ 
brane  in  this  region.  In  the  outlet  region  of  gases,  the  current 
densities  decrease  slightly  along  the  flow  channel  due  to  the 


2 

Current  density  (A/cm ) 

Fig.  3.  The  polarization  curves  for  various  relative  humidities  of  in¬ 
let  gases  in  the  coflow  (a)  and  counterflow  (b)  modes  (T  =  353  K, 
Pa  =  Pc  =  4  X  105  Pa,  £a  =  1.5,  £c  -  2.0,  aa  =  ac). 
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Fig.  4.  Distribution  of  water  content  of  membrane  with  dry  gases  in  the 
counterflow  mode  (T  -  353  K,  7av  —  1.0  A/cm2,  pa  =  pc  =  4  x  105  Pa, 
4  =  1.5,  $c  =2.0). 


Fig.  6.  Distribution  of  membrane  resistance  along  the  gas  channel  for  vari¬ 
ous  relative  humidities  of  inlet  gases  in  the  counterflow  mode  (T  =  353  K, 
4v  -  1.0  A/cm2,  pa  =  pc  =  4  x  105  Pa,  £a  =  1.5,  £c  =  2.0,  aa  =  ac). 


existing  of  liquid  water  in  the  coflow  mode.  As  discussed 
above,  liquid  water  lowers  the  effective  porosity  for  oxy¬ 
gen  gas  transport  and  increases  the  cathode  overpotential. 
Compared  to  the  coflow  mode,  counterflow  mode  improves 


the  current  density  distribution  with  dry  or  low  humidity 
gases. 

Fig.  6  shows  the  distribution  of  membrane  resistance 
along  the  gas  channel  for  various  relative  humidity  of  gases. 
Current  density  is  kept  constant  at  1.0  A/cm2.  The  results 
show  that  the  membrane  resistances  are  high  in  the  inlet 
and  outlet  regions  of  hydrogen  due  to  low  water  content  of 
membrane  with  dry  or  low  humidity  gases  in  counterflow 
mode. 

Fig.  7  shows  the  distribution  of  fraction  of  water  produced 
that  transports  through  the  anode  along  the  gas  channel  for 
various  relative  humidities  of  gases.  It  can  be  found  that  the 
net  flux  of  water  across  the  membrane  is  not  in  the  same 
direction  over  the  entire  area  of  the  membrane  in  the  fuel 
cell.  The  different  directions  of  the  flux  of  water  across 
the  membrane  at  different  regions  in  the  PEMFC  results  in 
an  internal  humidification  of  the  reactant  gases  in  counter¬ 
flow  operating  with  dry  gases.  On  the  dry  hydrogen  inlet 
region,  where  humidified  oxygen  exits  the  cell,  the  anode 
side  of  the  membrane  is  dryer  than  the  cathode  side,  thus 


Fig.  5.  Distribution  of  current  density  along  the  gas  channel  for  various 
relative  humidities  of  inlet  gases  in  the  coflow  (a)  and  counterflow  (b) 
modes  (T  =  353  K,  F^eil  =  0.7  V,  pa  =  pc  =  4  x  105  Pa,  £a  —  1.5, 
£c  =  2.0,  aa  =  ac). 
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Fig.  7.  Distribution  of/ along  the  gas  channel  for  various  relative  humidi¬ 
ties  of  inlet  gases  in  the  counterflow  mode  (T  =  353  K,  7av  =  1.0  A/cm2, 
Pa  =  Pc  =  4  X  105  Pa,  £a  =  1.5,  £c  -  2.0,  aa  -  ac). 
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establishing  a  relatively  steep  water  concentration  gradient 
across  the  membrane  (see  Fig.  4).  Water  by  electro-osmotic 
drag  from  the  anode  side  to  the  cathode  is  lower  than  the 
back-diffusion  water.  It  results  in  a  net  water  flux  from  cath¬ 
ode  to  anode  and  the  hydrogen  is  internally  humidified.  In 
the  outlet  region  of  hydrogen,  where  the  dry  oxygen  enters 
and  humidified  hydrogen  leaves  the  cell,  water  produced  by 
electrochemical  reaction  is  not  enough  to  humidified  dry 
oxygen,  the  cathode  side  of  the  membrane  is  dryer  than 
the  anode  side,  water  is  likely  to  be  transported  through 
the  membrane  from  the  anode  side  to  the  cathode  side  by 
electro-osmotic  drag  and  diffusion  [9].  The  dry  reactant 
gases  can  be  internally  well  humidified  and  maintain  high 
performance  when  PEMFC  is  operated  in  the  counterflow 
mode  without  external  humidification.  The  more  the  anode 
is  humidified,  the  more  water  will  migrate  from  anode  to 
cathode  in  both  coflow  and  counterflow  modes. 

3.3.  Effect  of  pressure  differential 

Net  water  transport  in  the  membrane  arises  from  three 
mechanisms:  electro-osmotic  drag,  diffusion  and  convection. 
At  low  current  density,  water  content  gradient  in  the  mem¬ 
brane  is  not  steep  and  water  content  of  membrane  on  the 
anode  side  is  high  enough  to  keep  the  membrane  proton  con¬ 
ductivity  well.  At  high  current  density,  water  transport  by 
electro-osmotic  drag  is  high  and  back  diffusion  water  by  the 
concentration  gradient  of  water  is  insufficient  to  replenish 
the  water  loss.  Membrane  on  the  anode  side  will  be  dehy¬ 
drated  and  membrane  resistance  increases  in  this  case.  Pres¬ 
sure  differential  can  be  used  to  force  water  transport  from 
cathode  to  anode. 

Figs.  8  and  9  show  the  distribution  of  membrane  resis¬ 
tance  and  fraction  of  water  produced  that  transports  through 
the  anode  in  the  counterflow  mode.  While  the  anode  pres¬ 
sure  is  kept  constant  at  4  x  105  Pa,  three  cathode  pressures, 
4  x  105  ,  5  x  105  and  6  x  105  Pa,  are  evaluated.  The  results 
show  that  higher  pressure  differential  results  in  lower  mem- 


Fig.  9.  Distribution  of  /  for  various  pressure  differentials  along  the 
gas  channel  in  the  counterflow  mode  (T  =  353  K,  7av  =  l.OA/cm2, 
pa  =  4  x  105  Pa,  £a  =  1.5,  £c  =  2.0,  aa  =ac  =  0.5). 

brane  resistance  over  the  whole  region  along  the  flow  chan¬ 
nel  and  higher  cell  performance.  A  higher  pressure  differ¬ 
ential  results  in  higher  value  of  f  that  is  to  say,  more  water 
discharge  from  anode.  Therefore,  a  higher  cathode  pressure 
can  be  used  to  reduce  the  anode  dehydration  problem  and 
improve  the  fuel  cell  performance.  However,  it  is  possible 
for  the  thin  Nation  membrane  to  be  damaged  if  high  pres¬ 
sure  differential  used  is  too  high.  In  addition,  high  pressure 
differential  may  result  in  seal  problem  in  fuel  cell,  especially 
in  fuel  cell  stacks. 

3.4.  Effect  of  cell  temperature 

Fig.  10  shows  the  polarization  curves  for  three  different 
cell  temperatures  in  the  coflow  mode.  The  results  show  that 
fuel  cell  performances  increase  with  increasing  temperature 
in  range  of  313-353  K.  Water  content  in  the  anode  side  of 
membrane  increases  and  membrane  resistance  decreases 
with  increasing  back  diffusion  coefficient  of  water  in 
membrane.  In  addition,  liquid  water  limits  the  transport  of 
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Fig.  8.  Distribution  of  membrane  resistance  for  various  pressure  dif¬ 
ferentials  along  the  gas  channel  in  the  counterflow  mode  (T  =  353  K, 
7av  —  l.OA/cm2,  pa  =  4  x  105  Pa,  £a  =  1.5,  £c  =  2.0,  aa  =  ac  =  0.5). 


Fig.  10.  The  polarization  curves  for  three  different  cell  temperatures  in  the 
coflow  mode  (pa  =  pc  =  4  x  105  Pa,  £a  —  1.5,  £c  =  2.0,  aa  =  ac  =  0.5). 
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Fig.  11.  The  polarization  curves  for  three  thicknesses  of  membrane  in 
the  coflow  mode  (T  -  353  K,  pa  =  pc  =  4  x  105  Pa,  £a  =  1.5,  £c  =  2.0, 
aa  =  ac  =  0.5). 

reactant  gases  and  covers  part  active  site  of  electrode,  when 
cell  temperature  increases,  part  of  liquid  water  evaporates 
and  cell  performance  increases.  However,  cell  temperature 
should  not  be  higher  than  373  K  because  reactant  gases  are 
diluted  by  water  vapor  and  lifetime  of  Nation  membrane 
decreases  at  high  temperature. 

3.5.  Effect  of  membrane  thickness 

Fig.  1 1  shows  the  polarization  curves  for  three  thicknesses 
of  membrane  in  the  coflow  modes.  Cell  performance  in¬ 
creases  with  decreasing  the  thickness  of  membrane  in  the 
both  modes,  due  to  the  higher  water  content  of  membrane 
and  shorter  distance  of  proton  across  the  membrane.  At  av¬ 
erage  current  density  7av  =  1.0  A/cm2,  the  cell  potential 
varies  from  0.718  V  (for  8m  =  0.5  x  10-4  m)  to  0.301  V  (for 
8m  =  1.25  x  10-4  m).  Fig.  12  shows  the  water  content  pro¬ 
files  of  membrane  at  v  =  0.06  m  in  the  coflow  mode.  Cur¬ 
rent  density  is  kept  constant  at  1.0  A/cm2.  The  simulation 
results  show  that  the  water  content  of  membrane  increases 


Fig.  12.  The  water  content  profiles  for  three  thicknesses  of  membrane 
at  v  =  0.06  m  in  the  coflow  mode  (T  =  353  K,  7av  =  1.0  A/cm2, 
Pa  =  Pc  =  4  x  105  Pa,  £a  =  1.5,  £c  =  2.0,  aa=ac  =  0.5). 


with  decreasing  the  thickness  of  membrane,  and  water  con¬ 
tent  gradient  across  the  membrane  is  high  for  thick  mem¬ 
brane.  It  seems  that  the  thinner  the  membrane  to  be  chosen 
as  electrolyte,  the  better  the  cell  performance.  However,  the 
difficulty  of  preparing  membrane  electrode  assemble  using 
thin  membrane  should  be  considered. 


4.  Conclusions 

A  two-dimensional  model  for  a  proton  exchange  mem¬ 
brane  fuel  cell  was  developed  in  which  water  transport  in 
the  membrane  by  electro-osmotic  drag,  back  diffusion  and 
pressure  differential  are  included.  Effect  of  liquid  water  on 
the  effective  porosity  for  gas  transport  in  the  cathode  and 
cathode  overpotential  were  considered.  The  model  is  used 
to  investigate  the  effect  of  flow  mode,  operation  conditions 
and  membrane  thickness  on  the  water  transport,  water  con¬ 
tent  of  membrane  and  performance  of  fuel  cell. 

The  simulations  show  that  the  water  content  of  membrane 
and  current  density  is  very  low  in  the  inlet  region  of  gases 
with  low  humidity  gases  in  coflow  mode.  The  dry  reactant 
gases  can  be  well  internally  humidified  and  maintain  high 
performance  when  PEMFC  is  operated  in  counterflow  mode 
without  external  humidification.  Compared  to  the  coflow 
mode,  counterflow  mode  improves  the  current  density  distri¬ 
bution  with  dry  and  low  humidity  gases.  Counterflow  mode 
does  not  show  any  advantageous  while  the  reactant  gases 
are  high  humidified  or  saturated. 

The  performance  of  cell  can  be  enhanced  by  increas¬ 
ing  the  temperature.  Pressure  differential  can  also  be  used 
to  improve  the  cell  performance.  The  results  show  that  a 
higher  cathode  pressure  can  be  used  to  reduce  the  anode 
dehydration  problem.  It  is  predicted  in  the  model  that  the 
water  content  of  membrane  and  cell  performance  increases 
with  decrease  in  the  thickness  of  the  membrane.  These  stud¬ 
ies  will  be  useful  for  the  design  and  operation  of  practical 
PEMFC  stacks. 

Calculated  polarization  curve  compares  very  well  with 
experimental  results  and  calculated  values  for  /av  are  agreed 
with  the  experimental  results. 
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